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Abstract 

We study the rate of convergence of some recursive procedures based on some "ex- 
act" or "approximate" Euler schemes which converge to the invariant measure of an 
ergodic SDE driven by a Levy process. The main interest of this work is to compare 
the rates induced by "exact" and "approximate" Euler schemes. In our main result, 
we show that replacing the small jumps by a Brownian component in the approximate 
case preserves the rate induced by the exact Euler scheme for a large class of Levy 
processes. 

Keywords: stochastic differential equation ; Levy process ; invariant distribution ; 
Euler scheme ; rate of convergence. 

1 Introduction 

In a recent paper (see |Pan05] ) . we investigated a family of several weighted empirical 
measures based on some Euler schemes with decreasing step in order to approximate 
recursively the invariant distribution u of an ergodic jump diffusion process X = (Xt)t>o 
solution to a SDE driven by a Levy process. More precisely, let (^fc)fc>i be such an 
Euler scheme with sequence of decreasing steps (7fc)fc>i and let {rjk)k>i be a sequence of 
nonnegative weights. We showed under some Lyapunov-type mean-reverting assumptions 
on the coefficients of the SDE and some light conditions on the steps and on the weights 
that, 

1 " - 

i^n(w,/) = \ \ y]??fc/(^fc^i(^)) °°' T^{f) a.s., (1) 

r/i + . . . + r?n ^ 

for a large class of functions / including bounded continuous functions (see Proposition [1] 
below, or |Pan05j for more general results, e.g., when u is not unique). We obtained this 
result for two types of Euler schemes: the "exact" Euler scheme that is built using the 
true increment of the Levy process and some "approximate" Euler schemes in which the 
Levy process increments are replaced by an approximation which can be simulated. 
The aim of this paper is to study the rate of a.s. weak convergence of (un) toward v for 
these schemes and to devise some variants of our schemes which speed up this rate. This 
problem has been first studied, for strongly mean-reverting Brownian diffusions, by Lam- 
berton and Pages ( |LaPa02] ) when r/^ = 7„, and by Lemaire ( |Lem06] ) for more general 
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weight sequences (see also |LaPa03j and |Lem05| ). In particular, Lemaire established in 
|Lem06) that considering some more general weights does not improve the rate obtained 
with r]n = 7n (although some choices may improve the "sharp" rate). Following this re- 
mark and in order to limit the technical difficulties, we will focus on the case r/„ = 7n- 
However, we will assume that (Xt) solution to the Levy driven SDE (see ([2])) is a weakly 
mean-reverting stochastic process, i.e. that (Xt) satisfies a weaker Lyapunov assumption 
than in the previously cited papers. 

As a first result, we show that the rate induced by the Exact Euler scheme (Scheme (E)) is 
the same as that obtained for Brownian diff'usions, provided the Levy process has moments 
up to order 4. In particular, the best rate is of order ns (see Theorems [U and [3]) . However, 
in practice, this "exact" scheme needs the increments of the jump component of the Levy 
process to be simulated in an exact way. This is not possible in general except in some par- 
ticular cases (stable processes, compound Poisson process, Gamma processes,. . . ). That is 
why we need to consider some approximate Euler schemes built with some approximations 
of the jump component, especially when the Levy process jumps infinitely often on any 
compact time interval. 

The canonical way to approximate the jump component is to truncate its small jumps 
(Scheme (P)). This amounts to replacing this jump component by a compensated com- 
pound Poisson process (CCPP). For this type of approximation, the smaller the truncation 
threshold is, the closer the law of the corresponding CCPP is to that of the true jump 
component, but conversely, the higher the intensity of its jumps is. So, there is a con- 
flict between the approximation of the jump component increments and the complexity 
of its simulation procedure (when there are too many jumps). The choice of the trunca- 
tion threshold is the result of a compromise between these constraints. It is time varying 
depending on the sequence (7„) and on the Levy measure. When the jump component 
has integrable variation, we show that it is possible to find a compromise which preserves 
the best rate of the exact Euler scheme. We mean that it is possible to construct a step 
sequence (7„) and a sequence of truncation thresholds such that on the one hand, the 
best rate induced by this type of approximation is of order n a (see Proposition [2]) and on 
the other hand the mean number of jumps at each time step remains uniformly bounded. 
This implies that the algorithm has a linear mean-complexity. Otherwise, this constraint 
of simulation slows down the best achievable rate. In particular, when the local behavior 
of the jump component is very irregular. Scheme (P) provides some very slow rates of 
convergence. 

We propose to overcome this problem by adapting a work by Asmussen and Rosinski 
( [AsRoOl] ) in which it is shown that when the truncation threshold tends to 0, the small 
jump component of a one-dimensional Levy process has asymptotically a Brownian be- 
havior. It can be extended to d-dimensional Levy processes (see Cohen and Rosinski, 
|CoRo05] ) . We then construct another Euler scheme (see Scheme (W)) by a wienerization 
of the small jumps. For this scheme, the compromise between the simulation and the 
approximation of the jump component is less constraining. Actually, we show that if the 
jump component has 3/2-integrable variation, it is possible to preserve the rate of ns and 
to respect the constraint of simulation. Furthermore, if vr is symmetric in a neighbor- 
hood of 0, the preceding assertion is valid without any conditions on the small jumps (see 
Theorem [2] and Proposition [2]) . 

Before outlining the structure of the paper, we list some notations: 

• The set M,^ ^ of matrices with d rows and / columns and real-valued entries will be en- 
dowed with the norm ||M|| := sup||2,|<i}. |M3;|/|x|. 

• For X G M*^ and k e N, x'^^ denotes the element of {W'-)^ defined by xf^ = 
Xi^Xi^ ...Xi„ for every ii, . . . G {1, . . . , /c}. 
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For every C^-function / : M'^ i— > M and x, y G M.'^, we adopt the following notation : 

ii,...,ifee{i,...,fc} '1 ^'^ 



If D^f is bounded, we set 

PVl|oo= sup sup I (X)|. 

n,...,ifcG{i,.--,fc} i^eK'* '^■^n • • • '^■^«fc 

• We say that V : ^ is an EQ-function (for Essentially Quadratic function) if V 
is a C^-function such that limy(x) = +cxd when — > +00, |Vy| < CVV and D^l/ is 
bounded. 

• We set r„ = Ylk=i 7fc! and for s > 0, r.[f^ = X]fc=i li''^- 



In Section [2l we introduce the framework and the algorithm, and we recall a result of 
convergence of the sequence of empirical measures established in |Pan05j . In Section [3l we 
state our main results about the rate of convergence induced by the exact and approximate 
Euler schemes when the Levy process has moments higher than 4. Sections H] (resp. [5]) 
are devoted to the proof of these results in the exact case (resp. approximate case). In 
Section [6l we state a partial extension of the main results when the Levy process has less 
moments. Finally, in Section [Tj we propose some numerical illustrations of our theoretical 
results. 



2 Setting and Background on convergence results 

For a Levy measure tt on M', we denote by (Hp) the following moment assumption 
(Hp) : / lyl'^P TT (dy) < +00 with p > 1. 

We recall that a Levy process with Levy measure vr is 2p-integrable (see e.g. [BaMiReOlj . 
Theorem 6.1). In |Pan05| . we studied the convergence to the invariant measure for every 
p > 0. Here, we only consider the p > 1 case because our main problem is to observe 
the impact of the approximation of the jump component which only depends on the small 
jumps. 

Throughout this paper, we denote by {Xt)t>o a solution to the following SDE 

dXt = b{Xt-^ )dt + a{Xt- )dWt + K{Xt- )dZt (2) 

where h -.W^ , a -.W^ ^Mdi a-nd k -.W^ ^ Mdi are continuous with sublinear growth, 
{Wt)t>o is a Z-dimensional Brownian motion and {Zt)t>o is a locally square-integrable 
purely discontinuous M'-valued Levy process independent of {Wt)t>o with Levy measure 
TT and characteristic function given for every t > by 

]E|gi<«,Zt>| ^ j ^i<u,y> _i_^^^^yy vr((iy))] . 

We recall that {Zt)t>o is a CCPP if and only if vr is a finite measure and that, otherwise, 
it can be constructed as a limit of CCPP: let (u„)„>i be a sequence of positive numbers 
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converging to 0. Let Dn = {\y\ > Un} and let ((^t,n)t>o)n>i denote the sequence of 
processes defined by 

Zt,n := J2 ^^-1{AZ,6D4 - 1 ! ynidy) Vt > 0. (3) 

0<s<i 

For every n > 1, {Zt^n)t>o is a CCPP with intensity A„ = 7r(L'„) and jump size distribution 
jin{dx) = 1 Dn ^(d\ ■ Furthermore, "-^+oo^ ^ ^2 locaUy uniformly, i.e. 

E{ sup |Zt-Zt,„|n ^^^±^0 vr>o. 

Q<t<T 

Discretization of the SDE. We introduce three Euler schemes. Scheme (E) is con- 
structed with the exact increments of the jump component and is called the exact Euler 
scheme. Schemes (P) and (W) are approximate Euler schemes. In scheme (P), we truncate 
the small jumps and in scheme (W), we refine the approximation by a wienerization of 
the small jumps. 

Let (7n)n>i be a decreasing sequence of positive numbers such that lim7„ = and such 
that r„ — +00. Let (?7„)„>i be a sequence of i.i.d. square integrable centered M^-valued 

random variables such that = Ii. Finally, let (Z„)„>i, (Z,^)„>i and (Z^)„>i be 
sequences of independent M'-valued random variables, independent of (C/„)n>i satisfying 
— jC, —p C, — C —p 

Zn = Z^n^ Z^=Z^^^n and Z„ = Z„ + y^Q„A„ Vn > 1, 

_ p 

where (A„)„>i is a sequence of i.i.d. random variables, independent of {Z^)n>i^Un)n>ii 
such that EAi = 0, Sa^ = and E{Af^} = 0, and {Qn) is a sequence of / x / matrices 
such that 

iQnQn)i,j = / yiVjiridy). 

We then denote by (Xn), (X^) and (X^), the Euler schemes recursively defined by 
Xo = Xq =Xq =xeR'^ and 

Xn+i =X„ + -Tn+iHXj + ,/^a{XjUn+i + (E) 

K+l = K+ ln+lh{K) + ^Ji;^(T{Xl)Un+l + /^(^n (P) 
Xn+l + In+lbiX^) + ,/j;;^a{X^ )Un+l + (W) 

We denote by (^n), {^n) {T^ ) the natural filtrations induced by [X^] and 

(X„ ) respectively. 

— p 

Remark 1. Note that Z^ can be simulated if both the intensity and the jump distribution 
of {Zt^n)t>o can be computed. Its simulation time depends on the number of jumps of 
{Zt,n)t on [0, 7„]. Its mean is 7r(Z)„)7„. In order to ensure the linear mean-complexity of 
the algorithm, we ask in practice these means to be bounded, i.e. 

sup7r(Z)n)7n < +00. (4) 

n>l 

In scheme (W), Q„ can be computed by the Choleski method as an upper triangular 
matrix if QnQn is definite. Otherwise, we can compute the principal square root of QnQn- 

The associated sequences of empirical measures are defined by 

where (77^) is a sequence of positive numbers such that i7„ = Ylk=i % "^+°°) +00. 
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Remark 2. As already mentioned, the rate of convergence will be only studied in the 
case ?7fc = 7fc- However, in the proof, we will intensively make use of convergence results 
for more general weighted empirical measures. That is why Proposition [1] is recalled in 
quite a general setting. 

Let us pass now to the Lyapunov mean-reverting assumption. Let a G (0, 1] be a parameter 
relative to the mean-reversion intensity. Let r > be a parameter relative to the growth of 
the noise coefficients a and k. In the sequel, we assume that there exists an E'Q-function 
V such that 

Assumption (Sa,r) : < CV Tr{aa*) + {{nf < CV with r < a. 

Assumption (Ra) : {VV, b) < fi - aV with a > and /? G M. 

The first deals with the growth control of the coefficients and the second is called the mean- 
reverting assumption. These two assumptions imply assumptions (Sa,p,q) and (Ra,p,q) 
introduced in |Pan05j . Hence, we derive the following result from |Pan05] : 

Proposition 1. Let a G (0,1], p > I and r G [0,a). Assume (Hp), (Ra) and (Sa,r)- 

Assume E{|C/i|^^} +E{|AipP} < +oo and {n]nMn) nonincreasing. 

(a) i. Then, 

supi>„(Fi+""^) < +00 a.s. (6) 

n>l 

Hence, the sequence (^'„)n>i is a.s. tight as soon as p/2 + a — 1 > 0. 

a. Moreover, if k{x) ' ' = o(|x|) and Tr(cr(T*) + < CV^^""^^, then every weak limit 
of{i'n) is an invariant probability for the SDE ([2]). In particular, if{Xt)t>o admits a unique 
invariant probability v, then for every continuous function f such that f = 0(1^2+""^)^ 
hm !>„(/) =K/). 

n— >oo _ _ 

Hi. Furthermore, E{yP(X„)} = 0(r„) and ifa = l, sup„>i E{yP(X„)} < +00. 

(b) The same result holds for {i'n )n>i CL'^d {i?^)n>i. 

Remark 3. For schemes (E) and (P), the above proposition is a direct consequence of 
Theorem 2 and Proposition 2 of |Pan05j . We did not study scheme (W) in |Pan05j but 
it is straightforward to show that the proposition holds true with a similar proof as that 
used for scheme (P). Note that when a = 1, n 1— > E{yP(X„)} is bounded whereas when 
a < 1, i.e. when the intensity of the mean-reverting is weak, one only has a control of 
its growth. This induces some technicalities but has no significant infiuence on the main 
results. 



3 Main results 

In this section, we suppose that E|ZtpP < +00 with p > 2 (see Section [6] for an extension 
to p G [1,2)). Let A denote the infinitesimal generator of (ATj). A is given for every 
C^-function / with bounded second derivative^] by, 

Afix) = (V/,6)(x) + ^Tria*Dya)ix) 

+ j {fix + Kix)y)-f{x)-{Vf{x),K{x)y)l{iyi<iy)TT{dy). 

^ Note that for such function, Af is well-defined since E\Ztf < +00. 
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We evaluate the rate of convergence on some test functions g such that g = Af + C where 
C is a nonnegative real number and / satisfies the following assumption: 

(CP) : (i) / G C^(M'^) and f{x) = 0{V{x)) as \x\ +00. 

{ii) For k = 2,3,4, f is a bounded and Lipschitz function. 

(m) |V/(x)p = 0{V^{x)) as |x| +00 with e e [0,p/2 + a - 1 - r). 

Since is invariant for the SDE ([2]), we know that v{Af) = (see e.g. |Pag01| ) and then, 
i/[g) = u{Af + C) = C. It follows that it suffices to evaluate the rate when C = 0. 

Remark 4. For a jump diffusion like ([5]), we are not able to characterize simply the set 
of functions g which can be represented as g = Af + C with / satisfying (C^ ). However, 
in the case of Brownian diffusions processes, some important works have been done in 
that direction. Actually, in [iPaVeOl] . |PaVe03j and [PaVe06j . Pardoux and Veretennikov 
show that in a Sobolev framework, existence and unicity hold for the Poisson equation 
g — u{g) = Af where A is the infinitesimal generator of a positive recurrent diffusion. 
Moreover, in |LaPa02j . Lamberton and Pages show that when the diffusion is an Ornstein- 
Uhlenbeck process, the above equation can be solved in C^(M'^). 

For this class of functions, the global structure of the rates of convergence is elucidated. 
For Scheme (E), our main result is Theorem [TJ We show that for every sequence (7^), 
there exists a sequence (pn) such that {pnV^{Af)) converges weakly: a fast-decreasing 
sequence (7^) (in a sense being precised in Theorem [TJa)) leads to a CLT and a slowly- 
decreasing sequence (7^) leads to a convergence in probability to a deterministic constant 
(see Theorem[H6)). The rate (p„) is maximal for a "critical" choice of (7,^) for which both 
types of convergence occur simultaneously. In particular, if 7„ = 7in~'' with C e (0, 1], 
the best rate holds for C = 5 (see "Particular Case"). In this case, p„ is of order ns. 
As concerns the approximate Euler schemes, our main results are Theorem [2] and Propo- 
sition [2j In the first one, we describe the structure of the rate induced by Scheme (P) 
and (W) as a function of (7^) and of («„). When (u„) decreases "sufficiently fast" in a 
sense depending on the choice of the scheme, on (7^) and on the Levy measure, the result 
induced by Scheme (E) remains valid for schemes (P) and (W). Otherwise, the approxi- 
mation of the jump component dictates a slower rate of convergence. 
Theorem [2] can not be directly applied in practice because it does not specify whether the 
fundamental condition of simulation ^ is compatible with the theoretical results. This is 
the purpose of Proposition [2] in which we give the best possible rates for schemes (P) and 
(W) under condition ([4]) as a function depending on the local behavior of the small jumps. 
In particular. Proposition [2] clarifies the impact of the wienerization of the small jumps 
announced in the introduction and shows that it makes possible to preserve the same rate 
of convergence of the exact Euler scheme for a wide class of Levy processes (for which the 
exact simulation of the increments is impossible). 

Let / G C^(R'^). We define by 

H^z, X, y) = f{z + K{x)y) - f{z) - {Vf{z), K{x)y) 

and z I— > {z,x,y) is denoted by Hlx,y. Our first result is the following: 

Theorem 1. Assume that E|ZfpP < +cxd with p > 2 and that ([2|) admits a unique 
invariant measure u. Let a G (0, 1] and r > such that (Ra) and (Sa,r) o,re satisfied and 
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p/2 + a — 1 > 2r. If moreover, E{?7®^} = 0, E{|C/ipP} < +00 and rjn = 7n for every 
n>\, then for every function / : M"' 1-^ M satisfying (C^ ), 

(a)If^^^=^^&[Q,+^),^nK{Af) ^ N(^m,a) 



+00, ^Vn^^f) 



V -L n 

w/iere o-j = f (k*V/p(x) + / [f{x + K{x)y) - f{x)f'Ti{dy))v{dx) and, 



bi{x) + (t>2{x) + 4>3{x))u{dx) with = ^D^f{x)b{xf 

M^) = I \D^f{x)-Mx); {o{x)ur^ + ^D^f{x){a{x)ur^^uAdu) 



and, ct)3{x) = ]- j 7r{dyi) j Tr{dy2)H^-^'y^{x,x,y2) 



+ I Adyi){{VHl,^y^{x)Mx))+ j ^uAdu)D\Hl,J{x){a{x)ur^). 



Particular Case. Assume that 7^ = 7in ^ with C, G (0, 1]. Then, 

'VlI^u^[Af) ^ AA(o,<72). ifC = i 
< ^^^AA(7m,<T2). if C G [1/3,1) 

^ ^J^nCpJA/) m if C < 1/3 



++00 



where 7 = if C G (1/3, 1) and 7 = if C = 1/3- On Figure[Tl one represents C 1-^ 

where denotes the exponent of the rate. One observes that max^/i(C) = /i(l/3) = 1/3. 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 



Figure 1: Rate of convergence for polynomial steps 



Remark 5. Theorem [T] shows that the rate is the same as that obtained for Brownian 
diffusions. In particular, when «; = 0, Theorem [T] extends the rate results of |LaPa02j 
and [Lem06j to the weakly mean-reverting diffusions (a < 1), whose convergence to the 
invariant measure has been studied in |LaPa03j . 

Note that the condition E{?7f ^} = is not necessary for the convergence of the empirical 
measures but plays a role in the rate. Without this condition, the best rate would be of 
order n*, obtained for (" = 1/2 (see |LaPa02] in the case of Brownian diffusions). 



Let us pass now to the main results for the approximate Euler schemes. Let {uk)k>i 
denote the sequence of truncation thresholds and set Pnjr = X]fc=i 7fc i|y|<uj. lul^'^^idy). 
For s £ {2,3,4}, we introduce a new assumption (Ag) which is relative to the impact 
of the jump component approximation as a function of the steps and of the truncation 
thresholds: 



(AJ : — - >a^G[0,+ooJ and — = > /J^ € [0, +ooJ. 



(s) 

Since s i— > 4 is a decreasing function, s i-^ Og and s i— > /3g both decrease. This can be 
interpreted as follows: the constraint on (uk) decreases with s. 

For s € {2, 3,4}, we also introduce another assumption on the Levy measure that will be 
necessary to transform some tightness results in some convergence in distribution results: 

. 9x r . /i|«|<M. • • • 

A^ : For every ii, . . . , G {1, . . . , /}, , , , , — converges in R. 

^ J|s/|<«Jy|^(^y) ^^=^1 

For instance, the above assumption is satisfied if 7r(dy) = '4j{y)\i{dy) where V satisfies: 
there exists a G [1,1 + 2) such that ^ Co G when y — > 0. 

Throughout this paper, we will say that the Levy measure vr is quasi- symmetric in a 
neighborhood of if /{|y|<M} V^^'^idy) = for u sufficiently small. In particular, this 
assertion holds if vr is symmetric in a neighborhood of 0. 

We will also say that a real-valued random variable X is quasi- suhgaussian if there exists 
m > and o" > such that for all M > 

P(|X| > M) < P(|y| +m > M) with y~AA(0,cj2). 

Theorem 2. Let a G (0, 1], r > 0, p < 2 such that the conditions of Theorem {1\ are 
satisfied. Assume that E{|AipP} < +oo and that {uk)k>i is decreases to 0. 

(a) i. Scheme (P): Assume that (Ag) holds with s = 2. 

• If ctg = or Ps = 0, then the conclusions of TheoremUl are still valid for [v^ ). 

• If dig G (0, +cxd] and (3s G (0,+cxd), then {—(ffi'n (^/))n>i is tight with quasi- suhgaussian 

011, TT 

limiting distributions. 

• If dig & (0, +c«] and = +oo, then (— (iyi^n (^/))n>i is tight with bounded limiting 
distributions. 

a. Scheme (W): Assume that (Ag) holds. 

Then, the conclusions of {a).i are valid for {v^ {Af))n>i with s = 3. Furthermore, if tt is 
quasi- symmetric in a neighborhood of and (A4) holds, the conclusions of {a).i are valid 
for {v^{Af))n>i with s = 4. 

(b) i. Scheme (P): Assume that (Ag) and (Ag) hold with s = 2. Then, 
-^u^{Af) M(m/as - ms,{&f/(3s)'^] if G (0, +00] and /3s G (0, +00) 

t->n,-K 

F p p 

^^{Af) — > m/ctg — nig if ctg G (0, +00] and Pg G (0,+oo), 



Pn,T7 

with \m2\ < rh2 = ('^/2)||-D^/||oo / ||K|p(x)i^((ix) and m and a'j like in TheoremUl 
a. Scheme (W): Assume that (A|) and (A|) hold. 

Then, the conclusions of {b).i are valid for {D^ {Af))n>\ with s = 3 and a real number 



satisfying Im^l < 7713 = {d^ /6)\\D^ f\\ao f \\K{x)\\'^u{dx). 

Furthermore, if 7: is quasi- symmetric in a neighborhood of and if {A\) and (A4) hold, the 
conclusions of {b).i are valid for {v^ {Af))n>i with s = 4 and a real number ru/j^ satisfying 
\m^\ <rh^ = (dV24)||L'Vl|oo / \\K{x)\\Mdx). 

Remark 6. Note that in the one-dimensional case, Assumption (Ag) is always satisfied 
when s = 2 or s = 4. In those cases, rUg = ^ f f^'^\x)K{xyh'{dx). If s = 3, Assumption 
(A^) is satisfied if /{|j^|<„^ y^TT{dy)/ /||y|<„^ \y\^7r{dy) ^ 03 G R. In this case, ms = 
aa^ / f^^\x)n{xYu{dx). In the multidimensional case, the value of is also explicit but 
its expression is more complicated (see proof of Lemma [T]) . 

Let us now state Proposition [2j In (a), we provide some conditions on the Levy 
measure in the neighborhood of which preserve the rate of convergence induced by the 
exact Euler scheme under the condition of simulation In (6), we suppose that the 
Levy measure has a density closed to that of an a-stable process in the neighborhood of 
and give in that case the optimal rate for the two schemes as a function of a. For these 
two parts, we also give some available choices of steps and truncation thresholds. 

Proposition 2. Let a G (0, 1], r > 0, p > 2 such that the conditions of Theorem\^ are 
satisfied. Assume that E{|Ai|2p} < +00 . 

(a) Assume that /||j^|<i} \y\'^'^{dy) < +00 with q G [0,2] and set 7^ = 71 A;^ 3 and Uk = Yk 
with r G Then, Condition ([4]) holds and, 

i. Scheme (P): If q < 1 and s = 2, nh!^(Af) ^ J2/3J\f(mV6, a'l) . 

a. Scheme (W): If q < 3/2 and s = 3, n^u^ (Af) ^ JlJ^N {m^/&,cr'^f) . 

Furthermore, if vr is quasi- symmetric in the neighborhood of 0, the preceding assertion is 
valid with s = 4 and every q G [0, 2]. 

(b) Assume that there exists eg > such that 

C C 

■K{dy) = ij{y)Xi{dy) with l|o<|y|<^^|^^ < -^{y) < (7) 

Set 7fc = 71A; 3 28~a' ^ Uk = 7^ with r [-, jjn^yj^]- Then, Condition ([4]) holds and, 
i. Scheme (P), s = 2: {n'^'^^hl{Af)) is tight, 
a. Scheme (W), s = 3; [rv^ <^-"'v^ {Af))^^^^ is tight. 

Remark 7. Figure [7| represents a 1— > h{a) where h{a) denotes the exponent of the 
optimal rate induced by each approximate scheme under the assumptions of Proposition 
Wip). This figure emphasizes the necessity of scheme (W) when the jump component has 
infinite variation because the optimal rate of convergence induced by scheme (P) decreases 
very rapidly in that case. 

Remark 8. The fact that we optimize the rate for the range of sequences (7fc, iifc)fc>i such 
that the linearity of the mean-complexity of the procedure is ensured can be disputable 
when the optimal rate is not of order ns . Actually, in this case, even if for a smaller level 
of truncation, the linearity of the complexity fails, the theoretical rate is better. Hence, 
another point of view consists in evaluating the order of precision as a function of the 
complexity. Some precise statements on that question would require some Berry-Esseen 
type estimates (in our inhomogeneous framework). Nevertheless, heuristic study can be 
done when ([7]) is satisfied and suggests that the asymptotic order of precision as a function 
of the mean-complexity is optimized for a class of steps and truncation levels including 
the choices of Proposition [2j 
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Figure 2: Optimal rate in terms of the local behavior of the Levy process 
4 Proof of Theorem [T] 

In this section, we prove the main result induced by the exact Euler scheme: Theorem [TJ 
Firstly, we decompose i>n{Af) (see Lemma [1|) and then, we compute the rate of each term 
of the decomposition in Lemmas [21 [3] and O (We will principally focus on Lemmas [2] and 
[3] where the rate of the jump part of the decomposition is studied). Finally, a synthesis of 
the previous lemmas is realized in subsection 14.21 and completes the proof of Theorem [TJ 

4.1 Decomposed computation of the rate of Un{Af). 
We set 

Xk,i = Xk-i + -fkK^k^i), and Xk^2 = -^fc,i + y/jkcr{Xk-i)Uk. 
Denote by (Z*'°')fc>i, a sequence of i.i.d. random variables such that Z*^' = Z and set 



Zh — Z, 



(fe) 

7fe 



fc=l 



Lemma 1. For f G C^(M'^), we have the following decomposition. 

n n 

^jkAf{Xk-i) = f{Xn) - f{Xo) - ^ (ei(7fc,X,.„i,?7fc) + e2(7fc,^fc-i,^'") 

k=l 

n 

- (Bi (7fc , 1 ) + 62 (7fc , Xk- i,Uk) + e3{jk,Xk-i,Xk,2,Z 

k=l 
n 

~Y,({Ri + R2){7k,Xk^i,Uk) + R3{lk,Xk.i,Xk,2,Z 



(k). 



k=l 



where, 

Uj,x,Uk) = ^{Vf{x),a{x)Uk), 

Ul,x,Z) = j\vf{x),K{x)dZs) + { Hf{x,x,AZs)-^ j Hf{x,x,yy{dy)), 



0<s<7 



©1(7, ^) = 7 / (v/(x + e^b{x)) - v/(x), h{x))de, 

Jo 

62(7, x,u)=j [ (1-9) {D^fix + jbix) + e^aix)u) - D^f{x)) ia{x)u)^^de, 
Jo 

10 



63(7, a;, Z) = Yl {H^ix,z,Zs-,AZs)-Hf{x,x,0,AZs)), 

0<s<7 

Ri{j,x,Uk) = ^{Vf{x + ^b{x)) - Vf{x),a{x)Uk), 
R2{l,x, Uk) = ^{D^f{x){a{x)Ukf' - E{D^f{x){a{x)Uk)^^}), 

Rsij,x,z,Z) = f\vf{z + K{x)Z,-)-Vf{x),K{x)dZs). 
Jo 

Proof. We write 

f{Xk) - f{Xk-i) = - f(Xk-i)) + (/(^l) - f{Xl)) + {f{Xk) - f{Xk,2) 

We expand the first two terms by the Taylor formula and use the Ito formula (with jumps) 
for the last one. The lemma follows by summing up the equality for k = 1, . . . ,n. □ 

As mentioned before, we study successively the rate of convergence of each term of the 
previous decomposition. We start by showing a CLT for the terms associated with ,^1 and 

Lemma 2. Assume that (Hp) holds for p > 2. Let f : ^ R satisfy (C^ ). Then, with 

the notations of Lemma [7], we have 

(a) 

EiMl, ^, Z)\^} =1 j {fix + <x)y) - f{x))\{dy) (8) 
and there exists (5 > and a locally hounded function C such that 

E{|6(7,x,Z)|2{i+^)}<C(x)7 (9) 
(b) Moreover, if (Ra) md (Sa,r) hold with 2r < p/2 + a — 1 and E{|[/ip^} < +00, then, 



1 " 

=y2Ui{ik,Xk-i,Uk)+uik,Xk-i,z^'^)) ^ M{o, 

^ k=l 



with a} = J (\a*Vf\^{x) + J {fix + K{x)y) - f{x))\{dy)) u{dx). 
Proof, (a). Let {Z, 

,n)n>i be the sequence of processes defined by ([3]). We know that 

E{ sup \Zs^n - Zs\'^} "^+°°) 0. 
{0<s<t} 

As Z^^n has bounded variations, ^2(7, 2:, -^.,11) can be written 

C2{'y,x,Z,^n) = ^ l{iAZ4>u„}{f(.x+K{x)AZs)-fix))--f (/(x+«;(x)?/) -/(x)) TT (dy) . 

o<s<7 J{\y\>'^n} 

Since D^f is bounded and E|Zf|^ < +00, one easily checks that ^2(7, a;, Z ^n) is a locally 
square-integrable purely discontinuous martingale. We deduce from the compensation 
formula that 

E{|6(7,^,^.,n)|'} = 7 / \f{x + K{x)y))-f{x)\Mdy). (10) 
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We also check that 



J {|y|<'f'i} 

Letting n — > +00 in (jlOp yields the first identity. 

Now, let us prove the inequality. ^2(7, 2;, Z) = {V f (x) , k(x) Z^) + where M is a 
martingale defined by 

M^= J2 [ i^fi^ + Ok{x)AZs) - V/(x), K{x)AZs)d9 
-if f\vf{x + eK{x)y)-Vf{x),K{x)y)de7r{dy). 







Let 6 G (0, 1] such that 4(1 + 6) < 2p. Since V/ is Lipschitz continuous, we derive from 
the Burkholder-Davis-Gundy inequality that 

E|Z^|2{i+^) <e{( |AZ,|2)^+'} and E\M^\^^'+'^ < C{x)e[{ ^ |AZ,|4)'+'}. 

0<s<7 0<s<7 

It follows that 

E{M7,x,Z)f <C^{x)E[{ J] \AZ,\'Y^'} + C2{x)e{{ J] |AZ/)^+^}. 

0<s<7 0<s<7 

Then, it suffices to prove that 

e|( I^^^P^'')^^'^} =0(7) for /> = and p = 2. (11) 

0<s<7 

Denote by (Ms) the martingale defined by Ms = Eo<5<7 lAZ^p+P - 7 J |y|2+P7r(dy). By 
the elementary inequality 

Vu,v G M and a > 0, |n + < 2"'^i-^(|u|" + (12) 

we have. 



e{( Y |AZ,|2+^)'+'} <C(E|M^|i+^ + 7i+^ /"|y|2+^7r(d2/)). 

0<s<7 

By the Burkholder-Davis-Gundy inequality, 

E\M^\'+^ < Ce[{ Y |AZ,|2(2+p))^|. 



i+s ■ 

AZs\''^'^P>)~ 

0<s<7 

Since (1 + 5)/2 < 1, it follows from (jl2|) and from the compensation formula that 



E|M^|1+^<Ce{ ^ |AZ,|(2+P)(i+^)| <C7 /"|y|(2+p)(i+^)7r(dy). 
Since 2 < (2 p){l + 6) < 2p, J |y|(2+p)(i+<5)vr((i?/) < +00. follows. 
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(6) Let {(Cfc )) k = 1, . . . ,n,n > 1} he a sequence of triangular arrays of square-integrable 
martingale increments defined by 



V ^ n 

Since S^/j = we have 

E{|a(7fc,^fc-i,f^fc)lV-^fc-i} = 7fck*V/|2(Xfc_i). 

Moreover, ^1(7^, C/fc) and ,^2(7/^ -^fc-i, -Z'*''') are independent conditionally to Tk-i 

and 

IE{ei(7fc,Xfc_i,C/fe)/^fc_i} = E{6(7fc,^fc-i,^*'V-^fe-i} = 0. 
Then, we deduce from ([8]) that 

E{l4lV-^fe-i} = T^(E{|ei(7fc,^fc-i,^:/fc)lV-^fc-i} +IE{|6(7fc,^fe-i,^*'')lV^^^ 

= ^(\a*Vf\HXk-i) + J {f{Xk-i + K{Xk-i)y) - f{Xk-i))\{dy)). 

Since f is bounded, we derive from Taylor's formula and from the assumptions on r 
and on V/ that 



(/(. + K{.)y) - f{.))\{dy) + |^*V/|2 < cy(^+-)^(2r) ^ (13) 
Hence, Proposition [1] yields 

f;E{|4"Pm.-i}^^^ / (k*V/|2+ / + -/(.)) V(dy))dz.. (14) 

J, 1 J 'J 



k=l 



Then, the lemma will follow from the central limit theorem for arrays of square-integrable 
martingale increments (see Hall and Heyde, [HaHeSOj ) provided the Lindeberg condition 
is fulfilled, i.e. 

n 

i?^ = ^E{|C^|2l^l^„l>^}/^,_i}^^^±^0 a.s. Vp>0. 

k=l 

Let A G (0, +00) and set 

n 

K,i = Yl '^{\x,-,\<A}^{\^kf'^mj^\>p}/^k~i}, 

k=l 

n 

Rn,2 =Y^{\X,_,\>A}^{\Ck\^M\i^\>p}/^k-l}- 
k=l 

We have E{\(^\H{^^n^^py/J^k-i} = F^iXk-ulk) where 

J- n 

Let 6 > such that ([9]) holds. By setting p = 1 + 5 and q = we derive from the 
Holder inequality that 

F^(x,7) < :f^E{|6(7,^,C^i)+6(7,^,^)l'^'+'^}^(niei(7,^,f^i)+6(7,^,^)l >/0V^) 

J- n ^ 
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On the one hand, we deduce from (I12p and from ([9]) that 

IE{|6(7, X, Ui) + 6(7, X, Z)|2(i+^)}TTi < C{x, 5){j'+^ + 7)TT7 < Ci{x, 

where x i— > Ci{x, 6) is locally bounded. On the other hand, we deduce from the Chebyschev 
inequality that, 

s -. ^ 
(F{\^i{j,x,Ui)+U-f,x,Z)\>py^))'^' < ^E{|ei (7, + 6(7,^,^)1'}^ 

<C2{x,6,p){^)Th 

^ n 

where x i— > C2{x, 6, p) is locally bounded. Then, for every ^ > and p > 0, 



Now, we observe R''^^. From ([l3]), we have : E{|^^|7jr^._i} < CVl^{Xk-i) with /5 < 
p/2 + a — 1. Therefore, 

<2 <El{|X, ,|>A}E{|4"lV-^.-l} < sup SUpP„(yi+"-^) = <A(A)^n(^^+"-^) 

where (/'(^) ^^^°°> 0. Since sup^gp^ < +oo (see Proposition [1]), letting 

A +00 yields 



□ 



Lemma 3. Let a e (0,1], r > and p > 2 suc/i that (Hp), (Ra), (Sa,r) /io/c? and 
p/2 + a - 1 > 2r. Assume that E{C/f ^} = and i/iai E{|C/i|2p} < +00. Let f -.W^ ^M. 
satisfying (C^ ). Then, 

(a) //ri'VVr^^o, 



1 

-= j;G3(7fc,Xfc_i,X,,2,^''') ^ 0. 

(h) //rf /Vr;^7e (o,+oo], 

1 " (fc) p /" 

-7rry]03(7fc,-^fc-l,-^fc,2,^ ) — ^ / </)3(x)z^(dx) 
J-n fc=l 

where (p3 is defined like in TheoremUi 



The proof of this lemma is realized in subsection | 

Remark 9. If Z is a compensated compound Poisson process, computing the rate of 
convergence of 63 consists in evaluating what happens after its first jump. Naturally, 
this argument has no sense when the Levy measure is not finite but the proof and the 
formulation of (j)^ show that it keeps some sense in average. 
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Lemma 4. Let a e (0,1], r > and p > 2 such that (Hp), (Ra), (Sa,r) hold and 

p/2 + a-l> 2r. Assume that E{Uf^} = and that E{\Ui\'^p} < +oo. Let f : R'^ ^ R 
satisfying (C^ ). Then, 

(a) ifv^^^ /Vr~n^o, 



^y]®i(^fc'^fc-i)+®2(7fc,^fc-i,^^fc) 0, 

n 

1 V — > — — — — (fc) P 

-=} Ri{-jk,Xk-i,Uk) + R2{'lk,Xk-.i,Uk)+R3{7k,Xk-i,Xk,2,Z ) — > 0. 



fc=l 



f6j //ri^^/-v/r^ ^ g (0,+oo], we have 



1 " 



.(2) 
n fc=l 



mi 



1 " 

and 02(7fc,^fc-i, ^fc) 



with mi = ^j D^f{x)b{x)^^i^{dx) 



and, m2=/ / ^I)V(^); K^;); (^^(x)^.)®^ ^ ^DV(x)(<T(x)nf4Pc/i (d7x)Kdx) 



1 " _ 

—72)" ^ -Rl(7fc) -'^fc-li t^fe) + -^2(7^:1 -'^fc-li ^^A;) + R3i'yk,Xk-.l,Xk^2, Z 



T.(2) 

^ " k=l 



0. 



Proof. The arguments of this proof are quite similar to those of the previous lemma. 
Then, we leave it to the reader. □ 



4.2 Synthesis and proof of Theorem [T] 

• Proof of Theorem [T] when Tn /^/Tn 0: Looking into the decomposition of Uni^f) 
introduced in lemma [H we deduce from lemmas [21 [3^a) and Hl^a) that 



(15) 



Now, / < CV. Then, by Proposition ^a).iii and Jensen's inequality, E{/(X„)} < 
E{VP{Xn)} < CT^. It implies that 

f{Xn) - f{Xo) Ll 



n— >+oo 







and Theorem [T] is obvious. 

• Proof of Theorem [1] when Tn /\/T\i — > 7 ^ (0, +00]: in this case, vTn 

< err. It 

implies that 

fjXn) - fjXo) 

p(2) n— >+( 

According to Lemmas [21 [3)^6) and[l|^6), we have 



0. 



- /-.N, ( f{Xn) - f{XQ, 



r(2) 



T.(2) 
J- n 



n— >+oo 



if 7 = +00 



n— >+oo 



M['ym,a'j^ if 7 < +00 



and the result follows. 
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4.3 Proof of Lemma [3] 

In the proof of Lemma El we usually need to show that some sequences tend to in 
probability. The arguments used for this are collected in the following lemma (these 
arguments also work for the proof of Lemma . 

Lemma 5. Let a £ (0,1], r > and p > 2. Assume (Hp), (Ra) o,nd (Sa,r)- Suppose 
that E{|C/ipP} < +CXD and let (Fi^) be a sequence of random variables such that is J-^- 
measurable. 

(a) Assume that F^^^/yT^ 0. 

^. //iFfcl <C772yf+--i(X,,_i), then, l/VT^ELi^fc-i ^ 0. 

n — >-\-oo 

a. If¥.{Fk/n-i] = and E{|Ffc|V^fc-i} < C {^lVP{Xk-i) + ^Iv"^ {Xk-i)) with e G 
[0,1) then, l|^/T-^T^k=lFk ^ 0. 

(b) Assume that Tn'* /^JT^ n^+oo^ ^ ^ (0,+cxd]. Then, 

n— ++00 

a. If¥.{Fk/Fk-i] = and E{|Ffe|V^fe-i} < C(7|y^P(^fc-i) + 7i^^(^fc-i)) with e e 
[0,1) then, iM'^ELi^fc ^ 0. 

Proof, (a) i. By Proposition [T]^a).iM, E{yP(X„)} < Cr„. We then derive from Jensen's 
inequality that 



^" fc = l ^''"k = l 



where p = p/2 + a — 1. Hence, the first assertion is obvious if 

If (fT6j) is not fulfilled, then we have liminf ^= X]fc=i7?\/r^ > because p/p < 1/2. It 

V-L n 



follows from the Kronecker Lemma that we have necessary X]A;>i7fc ~ +oo. By setting 
Vk = 7fc 1 we can apply Proposition [T] and deduce that 

1 " 

s^P^Y.^lV-^^''~\Xk-i)<+oo a.s. (17) 



Since /\/Tn 0, the first assertion follows when (|16p is not fulfilled. 

ii. Since E{yP(X„)} < r„, we derive from Jensen's inequality that E{|Ffcp} < C(7|rfc + 

7|r^), with e G [0,1). On the one hand, one checks that Yl^=ilk^k < (Ln )^- Hence, 
since Fn^ /^/Tn "'^'^°°> q, we have 



T^E^^'r,<ci^^^^0. (18) 



2r5 ^,2 



On the other hand, one observes that Ylk=i \ — Sfe=i — i^)" 2-e ^ +00. Hence, the 
Kronecker Lemma implies that 

^E^irf^^o. (19) 
fc=i 
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It follows that Y- X^fc=i l^ll-^feP} °°> 0. This yields the second assertion of (a). 
(b) i. We derive from the assumptions that 



n ^ n 

iyEl^^l^iyE^r^^^'^-nx.-O. (20) 



Then, {b).i. follows from p!7j) which is still valid because Ylk>i^k ~ 
a. It suffices to check that 

' J:e{|F,|^)^+oo. (21) 

n ) k=l 

With the same arguments as in {a).ii, one checks that 

7^E^Eil^^l'}^7^E^^'(r(^^r + 7^ witheE[0,l). (22) 

l-L n j A:=l " ) k=l \^ ri ) k=l 

On the one hand, we deduce from the Kronecker Lemma that 



1 " 

K'^n ) k=l 

On the other hand, for every e G [0, 1) 

l-L n ) k=l \^ n ) k=l 

because Ylk=i^k^k < (Tn^)'^- Hence, we derive from ([22]) that 

l-L n ; k=l 

which completes the proof. □ 

Proof of Lemma [S], (a) In order to alleviate the notations, we prove the lemma in the 
one-dimensional case. We set 

93(7, X, z, Z) = J ds J ■K{dyi){Hf [z + K{x)Zs,x,yi) - {x,x,yi 

©3(7, x,z,Z) is the compensator of 63(7, x, z, Z). Then, since 63(7, x,z,Z) is a purely 
discontinuous process, we have 

E{|(e3-e3)(7,x,z,Z)p} =E{^^(is j T:{dyi)\Hf {z + ^{x)Zs,x,yi) - Rf {x,x,yi)f}. 
By Taylor's formula, 

Hf{z, X, yi) = C {f'{z + eK{x)yi) - f{x)) K{x)yide. 
Jo 
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It follows that, 

1^-^(2; + K{x)Zs,x,yi)-H^ {x,x,yi)\ 

< sup \f{z + K,{x){Zs+6yi)) - f'{x + eK{x)yi)\.\K{x)yi\ 
0e[o,i] 

+ \nz + K{x)Zs)- nx)\.\K{x)yi\. 

f is a Lipschitz continuous function. Then, by setting z = x + 'yb{x) + ^/7(t(x)u, we 
deduce from the assumptions on the coefficients and from the fact that E{|Zsp} = 0(s) 
that 

E{|(e3-e3)(7,x,z,z)|2} 

<Cnfds j ^(dyi)(7'|6|2(x) + 7k(x)|V + |K(x)|2|Z,|2)|K(x)|V} 

< C(7V"+'^(x) +72(l + |u|2)y2r(^)). (23) 

Set Fk = (63 - e3)(7fc,^fc-i,^fc,2,^'''). Since ^{Fk/H-i} = 0, a + r < p and 2r < 
p/2 + a — 1, it follows from Lemma [5ll^a).ii and from the preceding inequality that 

1 " 

^(e3-e3)(7fc,^fc-i,^fe,2>^*'') ^0 when +00. (24) 

k=i 

Now, since /^^^ is bounded, we deduce from Taylor's formula that 

\Hf{z + K{x)Z,,x,yi)-Hf{x,x,yi)\<C\\K{x)f.\yi\\ 

Then, 

E{|e3(7,x,z,Z)|2} < C^^\\K{x)t j \yi\Mdyi? < Cj'V^^ix). 

By Lemma [5]^a).u, it follows that, 
1 

^E(®3(7fc,^fc-i,^fc,2,^''')-IE{e3(7fc,^fc-i,^fc,2,^''V-^fe-i}) ^^^0. (25) 
^ " fe=i 

Then, by (p4|) and (p5|) . (a) is obvious if we prove that 
1 " - 

-= VE{e3,i(7fc,^fc-i,^M'^''V-^fc-i} ^ and, (26) 

1 

— VE{e3,2(7fc,^fc-i,^fc,2>^''V-^fc-i} ^ with (27) 



©3,1(7,3;,^,-^) 

= 77f ds j %{dyi) dO (^f^^\z + k{x){Zs^ + %i)) - /(2)(z + 0K(x)yi)) (1 - e)K\x), 
63,2(7,^,^) =7 / ^(c^yi) f dO ^f(^){z + 9K{x)yi)-f^^Hx + eK{x)yi)) {1-9)k^{x) 



where 7f = f y1iT{dy){< +00) and vf is a probability measure defined by 7r((iyi) = yi-K{dyi)/Tt. 
Let us prove (p6]) . By Ito's formula, we have 

/(2)(z + K(x)(Z,^ + %))-/(2)(z + 0K(x)yi) = / /(3)(z + k(x)(Z,- +%i))K(x)dZ„ 

Jo 

+ ^ ^^*"(z + as:(x)(Z,- +0yi),x,AZ,). 

0<v<S7 

18 



Since /^^^ is bounded, the first term of the right-hand side is a martingale. Therefore, we 
obtain by the compensation formula that 

E{G3,i(7,^,^,^)} 

= -fTt [ ds U{dy{) [ de f^dv [ TT{dy2)E{Hf''\z + k{x){Z^ + eyi),x,y2)}il - e)K\x). 
Jo J Jo Jo J 

Finally, since 

\Hf'"' {z + K{x)iZ, + eyi),x, y2) \ < C\\f^^^ \\ooK\x)yl 

we have 

Since 2r < p/2 + a — 1, (j26|) follows from Lemma EJ^a).?. 

Now, let us prove ([27|) . Set z = x + ^b{x) + ^a{x)Ui. By Taylor's formula, there exist 

^1 G [x+^a{x)Ui+6K{x)yi,z+6K,(x)yi] and [x+9K{x)yi,x+^a{x)Ui+6K{x)yi] 
such that 

+ eK{x)yi) - + K{x)eK{x)yi) 

= lf^^Ki2)h{x) + ^f^^\x + eK{x)yi)a{x)U^+^f'^^\ii)a\x)Ul 

Since /^^^ and /^^^ are bounded and Ui is centered, 

|]E{G3,2(7fc,^fc-i,^fe,2)/^fc-i}| < C^l{\h\K\Xk-i) + a'^K^Xk-i)) < C^lV'^''''' {Xk-i) 

where we have used that a/2 + r < a V 2r. Since a V 2r < p/2 + a — 1, we deduce ([27|) 
from Lemma [5]^o).i. 

b) We keep the notations of (a). On the one hand, by ()23p and Lemma [5^6). m, we have 

n 

_ ^(63 - e3)(7fc,x,_i,Xfc,2, ^*") 0. 

^■n k=l 

On the other hand, we have 

E{|G3(7,x,z,Z)|2}<C7Vf(x) 
with e G [0, 1). Hence, by Lemma [5|^6).iz, it follows that 
1 " 

^ 5Z (®3(7fc, ^fc-i, Xfe,2, ^"") - E{e3(7fc, Xfc_i, Xfe,2, 0. 



Finally, it suffices to prove that 



^ VE{e3,i(7fc,^fe-i,^fc,2>^*")/-^fc-i} ^ "^3,1 (28) 



n fc=l 



1 " - 

and, Ve{ 93,2(7^,^^-1, ^fc,2,^*")/-^fc-i } ^ "^3,2. (29) 

n fc=l 
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with m^^i + m3^2 = / 4'3{x)i'{dx). In order to prove (|28p . we first show that 
1 " 

^^(E{e3,i(7fc,^fc-i,^fc,2>^'")/-^fe~i}-03,i(7fc,^fc-i,^fe-i,O)) ^ 0. (30) 



" A:=l 

Since /^"^^ is a bounded and Lipschitz continuous function, is also 2(5-Holder for every 
b e (0,1/2], i.e. 

[/^^]25= sup 1 <+00. 

It follows from the Taylor formula that 

^■^'^^ {z + k(x){Z^ + 6yi),x, y2)-H^^^^ {z + k(2;)(Z„ + 6*^1), x, 2/2) 

< c[/(^)]2.(|. - + n{xfW')<^?yl 
By setting z = x + 76(x) + ^a{x)u and taking (5 sufficiently small, we have 

W.{\Hf^'\z + K{x){Z, + eyi),x,y2)-Hf^'\x + K{x)9yi,x,y2)\] <C^\lM^^^ 
This implies that 



E 



{e3,i(7fe,^fc-i,^fc,2.^''') - e3,i(7fc,^fc-i,^fc-i,o)/.Ffc_i}| < c^l+W"2+-\Xk-i). 



Then, (j30p follows from Lemma [5][|6).i. 

Now, 03^i(7fc, Xfc„i, 0) is -measurable and 

|e3,i(7fc,^fc-i,^fc-i,o)| < C7fcV2'-(Xfc_i). 

Since 2r < p/2 + a — 1, we can apply Proposition [T] with r/^ = 7^. We obtain 

^ n -j^ " 2 /" 

-(^X]®3,i(7fci^fc-i,^fc-i,0) = -(^ X] "2^'^3,i(^fc-i) ""'"''°°> j (j)^^i{x)v{dx) a.s. 

with (/)3i(x) = ^ y 7r(dyi) j j T^{dy2)Hf''^\x + K{x)eyi,x,y2){l - 0)K^{x)yl 

_ 1 
~ 2 

It follows from (1301) that 



vr((iyi) / ir{dy2)H"-^'^^ix,x,y2). 



—(2j'^®sAlk,Xk-i,Xk,2,Z )/Tk-i] "~'+'^> m3,i = / <j)-i^i{x)v{dx) a.s. 

Finally, we prove ()29p . Since /''^•* and /*^^^ are bounded and Lipschitz continuous, we 
deduce from Taylor's formula that for every 6 S [0, 1/2], 

IE{63,2(7fc,-^fc-l,-^fc,2)/-^fc-l} = 7|(03,2(-^fc-l) +hfi{Xk-l)) + Pl{Xk-l,lk) + P2{Xk-l,lk) 

with 03,2(x)= / 7r(dyi) !\ef^''\x + eK{x)y{)b{x){l-e){K{x)yif 



Mx) = / 7T{dyi) / d9 / Pt;,(d^/)/(4)(x + 0/^(x)yi)^2(^)^2(^_^)(^(^)y^)2 



|pi(x,7)| < [/(')]n'+'|6(x)|'+'|K(x)|2 and |/^2(x,7)| < [/(')]257'+V(^)l'^'"*"'V(^)l'- 
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Since (a/2 + r) V (2r) < p/2 + a — 1, one can find 5 > such that 

|6(x)|^+^|k(x)|2 + |a(x)|2(i+^)|K(x)|2 < yf+'^-i(x). 

On the one hand, Lemma\^b) .i and the assumptions on the coefficients allow us to con- 
clude that 

1 ^ 

-(^^ (/Oi(Xfc_i,7fc) + p2(-^fc-i,7fc)) "^"^°°> 0- 

fc=l 

On the other hand, as |6|(x)|Kp(x) + |(Tp(x)|Kp(x) = o{V^~^'^^^ (x)) when |x| —>■ +oo, we 
derive from Proposition [1] applied with rjf^ = 7^ that 



n— ++00 

> 7713,2 



-^j^^liMXk-l) + 03,3(Xfe_i)) 
k=l 

with ?n,3,2 = J i4'3,2ix) + (p3,3{x))i'{dx). Checking that 

M^) + 4>3,3{^) = I TT{dyi){{Hl,^J{x)b{x) + I ¥u,{du){Hl,^y^f\x){cj{x)uy 
completes the proof. 



5 Proof of Theorem [2] 

The proof is built as follows. Like in the proof of Theorem[Tl we firstly decompose i>^{Af) 
and (Af) (see Lemma [6]). Some new terms appear due to the approximation of the 
jump component. That is why in the sequel, we focus on these parts of the decomposition 
(see Lemmas [3 and [8]). The other terms can be studied by the same process as their 
corresponding terms in the decomposition of z/„(^/) and then, are left to the reader. We 
denote by (z'*"'^ ),fc>i, a sequence of independent and cadlag processes such that 

{zl'^^^)t>o = {Zt,k)t>o and Z^^'^ = Vfc > 1 
For a C^-function / such that D^f is bounded, we define A^'^ and A^'^ by 

A'^Pfix) = (V/, 6)(x) + lTria*D^fa)ix) + [ {x, x, y)7r(dy) 

A^^'^fix) = A^'Pfix) + \ I D'fix)iKix)yrMdy), 

^ ■J{\y\<uk} 

These operators correspond respectively to the infinitesimal generators of 
dXt = b{Xt- )dt + a{Xt- )dWt + K{Xt- )dZt^k and 
dXt = b{Xt-)dt + a{Xt-)dWt + ^{Xt-)d{Zt,k + QkWt), 

where is a g-dimensional Brownian motion independent of W and (Z^ fc)t>o. 

Lemma 6. For a C"^ -function f such that D^f is bounded, we have the following decom- 
positions 

n n 

1) ^7fc^/(^r-i) = +/(^n) - /(^) - E {ii{ik,xi_,M+e2A^k,x',_,,z''''')) 

k=l k=l 
n 

- E (0i(7fc,^r-i) + ©2(7^,^^-1, + e3(7fc,^r-i,<2,^'''")) 

k=l 

n 

-Y,{{Ri + R2){ik,xl^iM + R3{ik,xl_^,xl^,z''"''')), 

k=l 
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where X^^ = xl_^ + T^K^I-i) + Vrk<K-i)Uk, = Ek=i^kiAf - A^^^ f){xl_^) 
and for a cddldg process Y , 

iip{l,x,Y) = f\vf{x),K{x)dYs)+{ V #/(x,x,Ay,))-7 / Hf{x,x,y)TT{dy)). 

n n 

2)^j,Afix:_,) = gZ + jZ + fixZ) - m - E (ei(7.,^ri,^.) +ct,B(7.,^ri,^'^'") 

k=l k=l 

n 

k=l 

n 

k=l 

where X^^^ = X^_^ + + 7fc6(X^__ J + ^a{X^_^)Uk, 

n n 

GZ=Y.ik{Af-A^^^f){x:_,i j: =-E/(^r)-/(C3)+7.(^'^'^/-^'='^/)(^ri), 



fc=l k=l 



with 3 = Xj^ - K(X^_i)QfcAfc. 

The two following lemmas are devoted to the additional terms of the preceding decompo- 
sition. In Lemma [71 we compute the rate of G„ and G„ and in Lemma [U we show that 

w 

does not have any consequences on the rate of the procedure. 

Lemma 7. Let a G (0,1], r > and p > 2 such that (Hp), (Ra), (Sa,r) hold and 
2r < p/2 + a - 1. Suppose that E{\Ui\'^p} < +oo and E{|Ai|2p} < +oo. Let f : R'^ ^ R 
satisfying (Cf). Then, 

(1) I. //lim„^+oo /3S < S^=i - A^^Pf){xl_,) ^ 0. 

(2) 

a. If lim„_+oo /3n,7r = +00, 

1 " -p 
limsup^E7fc|(^/-^''"^/)(^Li)| <^2 a.s. 

where fh2 = j \\^\\'^{x)v{dx). Furthermore, if {A\) holds, 

1 " 

-^Y.^k{Af - A'^'"' f){xl_^) m2 a.s. with jmsl < ma- (31) 



(2) Assume that s = 3 or that s = 4 if tt is quasi- symmetric in the neighborhood of 0. 

P ^ 



z. //lim„_+oo/3S<+oo^ELi7fc(^/-^'''^/)(^r-i) ^ 0- 



a. //lim„_,+oo/?i!]r = +00, 



1 " I 

hm snp ---Y,lk\{Af-A''''^f ) (X, 



w , 



< ms a.s. 



where = Cs||Z?''/||oo / ||'^(a^)||*i^(f^a^) with C3 = ^ and C4 = |j. Furthermore, if , 
(Ag) /loMs, 

^E^'^(^/-^'''^/)(^fc''-i)^^^^"^^ with|m,| <7fi,. (32) 
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Proof. By Taylor's formula, we have 

\Af{x)-A'''^f{x)\ = \f D'f{Q{^{x)yrMdy) < / |(^(x)y),(K(x)y),|vr(dy). 



For z eR'^,\ ZiZjl < \z\l < d\z\^. It follows that 

\Af{xl,) - A''^f{xl,)\ < ^\\D'f \\ooMxl,)f [ \y\Mdy)- (33) 

— P 

Since r < p/2 and E{yP(X^_^)} < Tfc, we deduce that 

n n „ 

^^,E{\Af{xl,) - A'^^f{r,_,)\} < Y,lk / \y\Mdy)Vn. 



fc=l k=l 



(2) 

Now, as lim„_»+oo /5n,Tr < +oo, the Kronecker Lemma yields 



^ n 

^E^W |y|Md2/)v^^-^o. 



The first assertion is obvious. 

a. Since lim„^+oo l^nlr = +oo, we deduce from Proposition[T]with r/^ = 7^ i|j^|<„^ lyP'^('^y) 
that 

4rE^^ / \y\Mdy)Mxl_,)f [ Mx)fi.{dx) a.s. 

,7r k=l •'\y\<'^k J 

because = 0(^2+"^^). Then, the second assertion follows from (I33p . 

Assume now that (Aj) holds. Since f is Lipschitz continuous, we deduce from Taylor's 

formula that 

Af{x) - A^^Pfix) = 

with pk{i,j)= yiyjTr{dy), ipi,j{x) =y^Ki^i-——Km,j{x)) 

JUvKu,.} r-' OlOm 



l,m 



^nd\Rl^{x)\<!^^y^^^^^\y\Mdy)\Hx)f. 

According to (A|), for every lim/)fc(i, |yp7r(dy) = € R. Set % 

■/'{|i/l<«fc} = Efc=i%- Then, 

Pn,iT k=l " 1 n 



k=l k=l 



with = {pkihj) — C(ijr]k)/r]f^. Firstly, since tj^ij < CV^ and r < p/2 + a — 1, Proposition 
□ applied with % = 7^. /{|y|<„j \y\'^7r{dy) yields 

1 " -p /■ 

j^'^Vki^ijiXk-i) 7pij{x)u{dx) a.s. 

" fc=i 

Secondly, = o(l). Since sup„>;^ l^fc=i Vky^^"'~'^{xl^_i) < +00 a.s., it is then easy 
to check that 



1 

jr^d^ki^i,ji^k-i) a.s. 
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The same argument is appropriate for i?f with el = (/||j^|<„^} \y\^T^{dy))/ (/{|j;|<„^} |yp7r((iy)). 
Finally, we obtain 

Pn,TT k=l i,j 

2) We derive from Taylor's formula 

\Af{x) - A'^'^fix)\ < Cs\\D'f\U\<xW I W<dy). (34) 



\y\<Uk 



with s = 3 and Cs = ^, or s = 4 and C4 = 54 if ijy|<„^ y®^7r((iy) = 0. As 2r < p/2 and 
E{yP(X^i)} < Tfc, it implies that 



I 1 ; 1 J\v\<Ui, 



k=l k=l 



\y\<Uk 

is) 



with s = 3 or s = 4 if J|j^|<„^ y'^^'^{dy) = 0. If lim„^+oo l3n,l < +00, we derive from the 
Kronecker Lemma that 

1 " r 

-^E^'^ / \y\'Ady)Vn^^^o. 

V J- n ^^-^ -'|s/|<«fe 

and the first assertion of (2) follows. 

Assume now that lim^^+oo /^n'^l = +00. Applying Proposition [1] to f{x) = \\k{x)\\^ with 
r]k = Ik ^\y\<^^ \yYT^{dy) yields 

1 V — ^ / —W 

limsup-7^2^7fe / \yY'K{dy)\\K,{Xj^_^)f < +00 a.s. 



Then, {2).ii follows from ([3^ 
Finally, the proof of ([52]) is similar to that of ([HT|) . □ 

Lemma 8. Lei a G (0,1], r > and p > 2 such that (Hp), (Ra), (Sa,r) hold and 
2r < p/2+a-l. Suppose thatE{\Ui\^P} < +00 and that E{\Ai\'^p} < +00. Let f ■.R'^^R 
satisfying (C^ ). Then, 

^=E(/(^r)-/(^M)+7.(^'='^/-^^'^/)(xr,)) ^ 0. (35) 

(2) — V / 

F/^ V F77, k=l 

Proof. Since 

(^A',Wf _ A^^Pf){x) = \f D^f{x){K{x)yr\{dy) = \e{D^ f {x){n{x)QkKf,r^] , 

we derive from Taylor's formula that 

fix:) - f{xi,)+ik{A'^'f - A'^^)f{x:_,) = 6(7.,<3,^r-i,Q.A,) 

~ , — w ,~ — w — w 

+ Rk,l{lk,X^__^, Qk^k) + Rk,2{'yk, X)^^^, Xf^_^,Qi:Ak) 
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where, ^1(7, z, x, QkAk) = ^{Vf{z), n{x)QkAk), 

Rk,i{^,x,QkAk) = ^{D^f{x){K{x)QkAk)''^ -nD^fi^){'^{x)QkAkf^}) 

i?fc,2(7,^,^,QfcAfc) = 7 {D^f{z + 9^K{x)QkAk) - D^f{x)){l - e){K{x)QkAk)'^^de . 

Jo 

Setting On = v r„ V , it suffices to show the three following steps: 
^)0-'El=iUlk,X^.„QkAk) ^ 0, 

h)e-^j:k=iRkAik,x^ii,QkAk) ^ 0. 
c) e-'Ek=iRkdik,x^,^,x^^i,QkAk) ^:^^o. 

a) We set 

^1(7, z, X, QkAk) = 6,1(7! X, QkAk) + 6,2(7, z, X, QkAk) 

with 6,1(7, a;, = ^{\/f{x),K{x)v) and 6,2(7, 2;, x, u) = ^/^{Vf{z) - V f {x) , k{x)v) . 
Let {Mn,i) and (M„_2) be the (^n)-niartingales defined by 

n n 

E— w ^ — > — w — w 

6,i(7fc,^fe--i,QfcAfc) and M„,2 = 2^ 6,2(7fc, -'^fc,3, -'^fc-i, Qfe^fc). 

k=l k=l 

We have to prove that e-^Mn,i "^^°°) and 0-^M„,2 0. 
According to (C^ ) and the assumptions on k, we check that 

< M >n,l . [ I |2 / , XT^£+a-l 



<^Y.^kl \y\\{dy)v'^^^-\x:_,). 

n. 7, 1 J\v\<ui, 



Now, by ©, sup„>i l/r„ELi7fc^^+'^~'(^ri) < +00 a-s- Since 4|<„^ ^ 

that < M >n,i /Ln "^^°°> a.s. Then, 



^ -|M„,i| < -^|M„,i| ^ 0. (36) 



Now, we turn to (M„^2)- Since V/ is Lipschitz continuous, it follows from the assumptions 
on the coefficients that 

mi,2{lk,xl„xt,.QkAk)?/Tk-i} <C [ \y\Mdy)(llV^'{Xk-i) + 7lV'hx^.i) 

with e < 1 and p = p/2 + a — l. Then, by a variant of Lemma [Ha).ii and {b).ii, one checks 
that 

b) -Rfc^i is very closed to R2 introduced in Lemma [Hand the arguments are similar. 

c) As D^f is bounded, one observes that 

E{\Rk,2{7k,X^,2,Xk-i,QkAk)\V^k-i} < C7feV'(Xr_i). 
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Then, since 2r < p/2, a variant of Lemmata) .ii and {b).ii yields 



1 " 



(37) 

By setting Ck,e{^) = x + ^k^{x) + ,y^cr{x)Uk + K{x){Zk + OQk^k), we decompose the 
integrand Rk,2 as follows: 

On the one hand, set = ^/JkQk-^k■ By Taylor's formula, we have 

{D'fix + jkKx)) - D^f{x))yf^ = D^fiiD-Hkbixy^yf 

where f{u);v;y'^'^ := J2i,j{'^D'^ fi,jiu),v)yiyj and £ [x,x + -fkb{x)]. Thus, since D'^f 
is bounded, we deduce that 

E{z?V(^r-i);7;^-K^r-i);('^(^r-i)Q;^Afc)^V-^ri} <Cjk I \y\Mdy)y^^^''{xti)- 

^ ^ J\y\<uk 

On the other hand, set /S.g{'^k-, x) = Ck,e{x) — (x + 'yb{x)). By Taylor's formula, 

(D^fiCkA^)) - D^fix + ^b{x)))yf = D^f{x)-^e{lk,x)-yf + ^l)V(ei); {^e{lk,x)r^-yf 
where f{u)-v'^^;y'^^ = T.^,j D^D^ fiAu))v'^^yiy, and H G [x + 76(2;), a,e(x)]. The 

— . W 

random variables t/^, Zj^ and are independent and independent of Then, since 

nUk/J'ti) = El^fc/^rj = n^f/^k-i} = 0' have 

E{z)3/(xri); A,(7,,xr_i); (K(xrii)QfcA,)^V-^ri} = 0. 
Now, since D^^f is bounded, one checks that 
E 



^ ^ Ay\<uk 

Since a/2 + r < p/2 + a — 1 and 2r < p/2 + a — 1, it follows that 

\E{Rk,2i7k,X^,2,xti,QkAk)/J'k-i}\ <C [ \y\\idy)jlV-2+'^-Hxti) 

By a variant of Lemma El^a).^ and {b).i, we derive from the previous inequality that 



n 

1 ^ — > ~ — W — W IP 

-— ^^E{i?,,2(7fc,X;,,2,^fc-i,QfcAfc)/^fc_i} ^ 0, 

V-L n V i n k=l 

Then, assertion c) follows from ()37p . □ 
6 An additional result 

In this section, we present a partial extension when the Levy process has a moment of 
order 2p with p £ (1,2]. In this case, stating some global results as in Theorems [1] and [2] 
would need two kinds of restrictions: either to assume that at least the derivatives of / 
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tend to when — > +cxd or to impose more constraints on the growth of the coefficients. 
The first alternative leads to a very technical proof and the second one can not be really 
envisaged for the drift term. Actually, we recall that in this type of problem, b produces 
the mean-reverting effect and then, it would not be natural to suppose that for instance, 
b is bounded . 

That is why we propose to state a partial result for fast-decreasing steps for which the 
extension does only require some weak restrictions on /. We introduce a new assumption 
on the steps depending on the intensity of the mean-reverting: 



t(2) 1 " aV(2r) 

nn^+oo^^ if„ = l and, -^J2jIT,^ n^+oo^Q ^ ^ ^_ ^3^^ 



Then, 



Theorem 3. Assume thatK{\Zt\'^P} < +00 with p £ (1,2] and that ([2]) admits a unique 
invariant measure v . Let a G (0, 1] and r > such that (Ra) and (Sa,r) are satisfied and 
such that p/2 + a — 1 > r. Suppose that E,{Uf^} = 0, E{|C/i|^} < +00 and r]n = 7n for 
every n > 1. Let / : R"^ 1-^ M 6e a -function having bounded derivatives and satisfying 
f{x) = 0{^V{x)) as \x\ +00. Then, 

(a) Scheme (E): If ^ holds, ^Vn{Af) ^ M(q,&'] 

(h) Scheme (P): If (p8]) holds and Pni-/\/Tn 0, the conclusion of (a) is valid for Scheme 
(P)- 

(c) Scheme (W): If ([38]) holds and PnJr/^/Tn — ^ 0, with s = 3 or s = A if n is quasi- 
symmetric in the neighborhood of 0, the conclusion of (a) is valid for Scheme (W). 

Remark 10. We refer to |Pan06j for a proof of this result. 

Assumption (j38p is less constraining when a = 1 because the L^'-control of the Euler 
scheme is better in this case (see Proposition [1]). Note that when a = 1, Theorem [3K a) 
corresponds to Theorem [T]( a) when 7 = 0. 

Let 7fc = 7iA;~^ with C, G (0,1] and 71 > 0. For Scheme (E), Theorem [3] applies in the 
following cases: 

C > - if a = 1 and C > — if a < 1 

^3 ^ 3p + 2?? 

where r\ = a\J (2r). Since \/r^ '^'k?^ j ^i^n-^ if C £ (0) 1)) we derive that for every 



-C 

e > 0, there exists an Euler scheme with polynomial step such that the rate of convergence 

1 _ V ^ 

is of order ns if a = 1 and n^p+a^? if a < 1. 



7 Numerical comparison of Schemes (P) and (W) 

1. When z^(/) can be theoretically computed. In this first example, we are interested 
in the two-dimensional SDE 

dXt = -Xt^dt + dZt (39) 

where Z is a symmetric purely discontinuous Levy process (having no drift term). We 
consider (j) : x |xp and denote by 1^, the unique invariant SDE ()39p . We can easily 
compute i^(</)). In fact, as vr is symmetric and i^iAcp) = 0, 

A^ = -2(l) + J \y\\idy) =^ '^{^) = ^ J |y|^vr(dy). 
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Let us test this theoretical result on a 2-dimensional example. Assume that 

7r(")(d2/) = l{\y\<i} ^2(rfj/) + ^{\y\>i}J^^2{dy) with a G (1,3). 

We have = 7r(l/(2 — a) + 1/4). In figures [3] and [71 we observe the rate for two values 
of Q taking the choices of steps and truncation thresholds of Proposition ^b). These 




Figure 3: n i-^ ^n{<P), a = 1 Figure A: n ^ Vn{<P), ol = 5/3 

simulations are coherent with the theoretical results. Indeed, when a = 1, the optimal 
asymptotic rates induced by Schemes (P) and (W) are the same (with order ns). When 
a = 5/3, the optimal asymptotic rate induced by scheme (W) is still of order na whereas 
that of scheme (P) is of order n^. 

2. Another example. Now, we observe the following two-dimensional SDE 

dXt = dt + (1 + \Xt-\)\dZt 

v 1 + \Xt- I 

where [Zt) is a purely discontinuous Levy process having no drift term with Levy measure 
vr^") (defined in the preceding example). One checks that Proposition [2] applies with 
V{x) = 1 + a = 3/4, r = 1/4 and every p S (2,3). As in the preceding example, 
we test our procedure in the cases a = 1 and a = 5/3. As the dynamical system is less 
stable, the convergence is slower but we can observe the same phenomenon. 




Figure 5: n i-^ ^n{(p), a = 1 Figure 6: n i-^ i^n(</'), a = 5/3 
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